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ABSTRACT 

We report on a BVRcIc survey of field W Ursae Majoris binary stars and present accurate 
colors for 606 systems that have been observed on at least three photometric nights from a robotic 
observatory in southern Arizona. Comparison with earlier photometry for a subset of the systems 
shows good agreement. We investigate two independent methods of determining the interstellar 
reddening, although both have limitations that can render them less effective than desired. A 
subset of 101 systems shows good agreement between the two reddening methods. 



Subject headings: binaries: close — binaries: eclipsing 



1. Introduction 

Overcontact binary starifl are systems in which 
both components share a common envelope be- 
tween the inner and outer contact surfaces, and 
are therefore in physical contact with one another, 
leading to luminosity and mass exchange. Such 
systems show light curves with very similar eclipse 
depths and little color variation over the orbital 
cycle. The W Ursae Majoris systems are the most 
well-known overcontact systems, consisting of F, 
G and K-type stars with convective common en- 
velopes. But early-type sy stems with radiative 
common envelopes do exist (jTerrell. et al. 1 12003) 



1 We note that there is a preference among some in the com- 
munity: theterm " contact" rathe r than "overcontact." 

iRucinski I dl997al) and I Wilson! J200lh discuss both sides of 
the issue. 



and a re even quite commo n among the O-type 
stars ijHilditch fe Bell Ill987h . 

Despite over four decades of efforts to un- 
derstand their internal structure and evolution, 
no model has emerged as clearly superior in ex- 
plaining th eir observ e d prop ertie s, althoug h activ e 
work (e.g., IStepieh I (|2009l ) and IStepieh I (|201ll) ) 
does show promise. Almost all of the theoreti- 
cal work on overcontact systems has involved the 
W U Ma systems, although the ev olutionary mod- 
els of iNelson fc: Eggletonl (|200lh show that the 
overcontact state is achieved for a wide variety 

of binary initial condi tions ( viz. their Figure 9). 

The seminal work by iLucv I (|l968l ) showed that 
overcontact systems with common convective en- 
velopes could exist at zero-age but the agreement 
of hi s models with the observed period-color rela- 
tion (Eggen 1967) was not very good and Eggen's 
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period-color relation has played a critical role in 
the testing of models of W UMa systems ever 
since. The primary goal of this survey is to pro- 
vide a significantly larger sample of overcontact 
systems with accurate photometry and to extend 
the photometry to the RqIc passbands. 

Although the proliferation of relatively inex- 
pensive CCD cameras and telescopes has resulted 
in a large body of light curve data for W UMa sys- 
tems, most of these data are instrumental differ- 
ential magnitudes rather than transformed stan- 
dard magnitudes, a circumstance dictated in most 
cases by the predominantly non-photometric skies 
available to the observers. iRucinski I (|2006h dis- 
cusses the period-color-luminosity relationship for 
W UMa systems and its use as a distance indicator 
for these systems. One goal of our program is to 
provide more accurate color data than is currently 
available for this calibration. Accurate colors are 
also necessary for light curve solutions where the 
spectral type of the system is not available. 

2. Sonoita Research Observatory 
2.1. Hardware 

The goal of our survey program is to measure 
accurate colors of a large number of W UMa sys- 
tems brighter than 14 th magnitude visible from 
our observatory in Sonoita, Arizona. The equip- 
ment at the Sonoita Research Observatory (SRO) 
consists of a 0.35m Schmidt-Cassegrain telescope 
on a robotic German equatorial mount, a Santa 
Barbara Instrument Group STL-1001E CCD cam- 
era and a Custom Scie ntific BVRcIc filter set 
with the iBessell I (|l990h prescriptions. With this 
equipment, we can achieve photometric precisions 
of 0.01 magnitudes at 14 th magnitude with rea- 
sonable integration times. 

The observatory sits on a private site located 
about 60 miles east-southeast of Kitt Peak Na- 
tional Observatory in very dark skies. Previously 
unknown main- belt asteroids in the 19 t/l -20*' 1 mag- 
nitude range have been discovered in five-minute 
unfiltered images taken at SRO, and stacked, unfil- 
tered images have revealed 23 rd magnitude stars. 
Seeing is typically 2"-3" full-width at half maxi- 
mum. 



2.2. Software 

As with all robotic observatories, the main de- 
terminant of observing efficiency at SRO has been 
the software used to control the telescope and 
camera. All of the observing software we use 
is commercially available. The two key software 
programs are ACP Observatory Control Software 
(hereafter, ACP) and ACP Scheduler from DC-3 
Dreams, Incorporated of Mesa, ArizonaH Over- 
all operations control is handled by ACP. At the 
beginning of each night, it opens the dome, and 
takes twilight flats. During the night, it monitors 
the weather conditions (closing the dome if the 
weather degrades to a dangerous state) , moves the 
telescope and dome, instructs the camera control 
software to take images, and monitors the tele- 
scope focus, refocusing as needed. At the end of 
the night, it takes more twilight flats and then 
parks the telescope and closes the dome. The en- 
tire operation is run unattended, allowing for a 
tremendous increase in observing efficiency over 
traditional human-driven telescope operations. 

The primary enabler of a survey project like 
this one is the scheduling software. The ACP 
Scheduler is a real-time scheduler that can ad- 
just to changing observing conditions and observ- 
ing priorities, as opposed to the traditional queue 
schedulers that decide the observing plan at the 
beginning of the night. Executing a survey like 
this one involves merely entering the coordinates 
of the desired fields (i.e. standard star fields and 
the object fields), filters and exposure times into 
the ACP Scheduler database once. At the start 
of each night, the scheduler chooses the best field 
to observe. When that field is completed, a new 
field is chosen based on the current conditions. 
If the weather degrades from photometric condi- 
tions, but is still reasonable for differential pho- 
tometry projects, the scheduler will begin doing 
those. Should the weather later improve, it will 
return to projects like the survey that require pho- 
tometric skies. The scheduler thus makes the most 
efficient use of the observing time, all without hu- 
man intervention. 



2 Web site at http://www.dc3.com 
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3. Data Gathering and Processing 

While the observatory operation for the sur- 
vey is untraditional, the data gathering is not. 
The goal is to measure the program stars in four 
Johnson-Cousins filters (BVRcIc) an d trans- 
form the observations to the standard system via 
the observation of Landolt (1983, 1992) standard 
fields. Our approach is to separate the determina- 
tion of the extinction coefficients from the trans- 
formation coefficients. We observe a reasonably 
bright star with a B — V value near zero through 
each filter periodically during the night as it moves 
through a reasonable range of airmass, typically 
from 1.0 to 2.5 airmasses. The star is chosen so 
that it can be observed low in the eastern sky at 
the beginning of the night, up to the zenith, and 
then in the western sky late in the night. This 
approach typically yields about a dozen measure- 
ments through each filter per night, and covers 
both sides of the meridian, providing a diagnostic 
of the photometric quality of the night. A linear 
fit of instrumental magnitude versus airmass then 
provides the (first-order) extinction coefficient. 

Observations of standard and target fields 
are made in the sequence Ic,Rc,V,B, and then 
B,V,Rc,Ic to minimize extinction issues. The 
two magnitude measurements for each star in 
each filter are averaged to provide coeval values 
for forming colors. Determination of the trans- 
formation coefficients is done in the usual way 
by observing several Landolt standard fields near 
the meridian during the night and then perform- 
ing linear fits to the various instrumental versus 
standard magnitude plots. 

As images are collected during the night, they 
are losslessly compressed and transferred to a stor- 
age server at the lead author's institution. Pro- 
cessing of the images is performed with scripts 
written in PyRAF, a Python interface to IRAF 
available from the Space Science Institute web 
siteH They are then bias and dark subtracted 
and flatfielded. Extinction coefficients are first de- 
termined and then the transformation coefficients 
determined. This part of the process is interac- 
tive, with plots that allow the user to determine 
the quality of the night and delete obviously bad 
images from the reduction. Once the extinction 



and transformation coefficients are known, the im- 
ages of target fields can be reduced and the ob- 
servations transformed onto the standard system. 
Aperture photometry is performed with an aper- 
ture radius of eight pixels, sky annulus radius of 
twelve pixels and an annulus thickness of three 
pixels. The pixels subtend 1".24 on the sky. 

4. Survey Targets and Results 

The sur vey began in 2005 with targets s electe d 



3 PyRAF is available at http : //www. stsci . edu/resources/sof tware. 



from the iPribulla. Kreiner. fc Tremko I (|2003l ). 
hereafter PKT, catalo g of 361 overcontact sys- 
tems. Soon thereafter, Gettel. Geske. &: McKavl 
(200G), hereafter GGM, published a catalog of 
1022 overcontact systems identified from ROTSE- 
I sky patrol images. The final list of observed tar- 
gets came primarily from these two catalogs, with 
two systems (GSC 0067-00348 and GSC 2545- 
00970) added after we made follow-up observa- 
tions of new discoveries. The last of the observa- 
tions reported herein were made in 2008. Observa- 
tions were made on 237 photometric nights, with 
this project sharing telescope time with several 
other survey projects. 

In total, 606 systems with 9 < V max < 13 had 
observations of sufficient quality to be listed in Ta- 
ble 1. The table is given in abbreviated form in 
the printed version of the paper to show the file 
structure, with the full table given in the online 
version. The star identifications are either the 
GCVS designation, if known, or the designation 
given in GGM. The objects listed by their inden- 
tification from the GGM catalog can be searched 
in the SIMBAD database by using the "ggm2006" 
catalog identifier. We chose to retain the GGM 
designation because SIMBAD does not list Guide 
Star Catalog identifications of many of these sys- 
tems. The positions, orbital periods and maxi- 
mum/minimum V magnitudes are taken from ei- 
ther the PKT or the GGM catalogs. In a handful 
of cases, the positions are from astrometry of our 
images where the GGM positions were ambiguous 
due to blending. All systems were observed on at 
least three nights but 251 systems had observa- 
tions on four or more nights. The B ~ V, V — Re, 
and Rc — Ic colors are given as measured, and the 
V — Ic color is formed as the sum of V — Rc and 
Rc ~ Ic- The errors given are the standard de- 
h 1 S*8titfi§'$fe-ific individual measurements, and this 
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include both the measurement errors and any in- 
trinsic color changes in the stars, which we assume 
to be small. 

To assess the accuracy of the colors, they were 
compared to values found in the Tycho catalog 
and from other sources in the literature. The 
bulk of the comparison data came from the Ty- 
cho catalog and were selected for comparison to 
our data if the error listed in the Tycho catalog 
was 0.05 mag or less. To extend the comparison 
past B — V = 1 , values for CC C om (jBradstreet I 
19851). VZ Psc dBradstreet Ill985l). an d V523 Cas 



(|Samec. Faulkner, fc Williams I 2004) were taken 
from the literature. Figure Q] shows our B — V val- 
ues plotted against the comparison data. Figure 
[5] shows our V — Ic values plotted against val- 
ues from Tycho and the literature. Once again, 
the red end of the figure was extended past the 
range of accurate Tycho values, wit h the value for 
XY Leo taken from Hrivnak I (| 19851) and values for 
CC Com, VZ Psc, V523 Cas, V574 Lyr, CT Tau, 
LR Cam, 766402, 7710246, and 13149917 taken 
from The Amateur S ky Survey Mark IV (TASS) 
(jDroege. et al. 1 120061 ). On the blue end, values 
for CT Tau and 2962451 were also taken from the 
TASS catalog. Both figures show good agreement 
between our photometry and that of others. 

We examined our data for any systematic de- 
pendences of the errors on various quantities. One 
obvious dependence would be on the magnitude 
of the binary. In a constant-exposure survey, this 
type of dependence is expected and easily revealed. 
Our survey, however, had exposures optimized for 
each individual system, and our targets are all 
in the brightness range where we expect intrin- 
sic color variability and the transformation to the 
standard system to be the dominant noise sources, 
rather than sources like photon noise or sky back- 
ground noise. In that case, we do not expect to 
see any strong dependence on the brightness of the 
system and Figure [3] shows this to be the case. 

Another possibility is a dependence of the color 
error on the color itself. Our exposure times were 
determined using the V magnitude of the binary 
at maximum brightness and the known sensitivity 
of the equipment for each filter. Redder systems 
would require a longer exposure in B than a bluer 
system of the same brightness in order to have 
the same signal-to-noise ratio for photon noise (the 
dominant noise source in a given exposure for our 




Comparison B-V 



Fig. 1. — SRO B — V values compared to B — V 
values in the literature for 35 systems. 




Literature V-l 



Fig. 2. — SRO V — Ic values compared to V — Ic 
values in the literature for 23 systems. 
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equipment and the magnitude range of our sur- 
vey). Since our exposures were set to achieve a 
signal-to-noise ratio of 200 or better (for photon 
noise) for each exposure, we still expect the ob- 
served errors to be dominated by intrinsic vari- 
ablity and transformation uncertainties. Figure [4] 
shows that this is indeed the case. 



5. Reddening 

To be of the most usefulness, the colors of the 
stars need to be corrected for the effect of inter- 
stellar reddening. Normally, an accurate mea- 
surement of the reddening requires a detailed 
spectroscopic a nalys is, such as that employed by 



Terrell, et al. I (I2003f) using the Na D lines as cal 



ibrated by iMunari fe Zwitter I (|19971) . For the 



majority of the systems in our survey, such data 
do not exist. We therefore present a preliminary 
analysis using two approaches for which complete 
data are available: determining a reddening-free 
quantity from our photometry, and computing 
the reddening from the neutral hydrogen column 
density, Nhi, along the line of sight to the sys- 
tem. The former approach, referred to hereafter 
as the color difference method, follows that used 
for early-type stars with U — B and B — V col- 
ors to determine the redden ing-free parameter 
Q ([Johnson fc Morga n, 1953T). The latter is the 



method employed by lRucinski I (jl997l ) in calibrat- 
ing the absolute magnitudes of W UMa systems. 
While these approaches have limitations, as dis- 
cussed below, we believe they can provide useful 
information until better data can be obtained. 

Following conventional notation, we write the 
observed color as B — V , the unreddened color as 
(B — V)q, and the reddening as E(B — V) so that 
B — V = (B — V)o + E(B - V). The color dif- 
ference method uses the difference in two colors 
to form a reddening-free quantity usually respre- 
sented by the letter Q. The key is to mutiply one 
of the colors by a constant factor representing the 
ratio of the color excesses for the two colors, i.e., 
Q = Cl- { -§^ x C*2. While the traditional Q 



quantity, (U — B) 



E{U-B) 



x (B — V), is not use- 



E(B-V) 

ful for the stars of spectral types in this survey, 
the basic idea can be applied to more appropriate 
filters. Here we use the B — V and V — Ic colors, 
and thus define Q as (B — V) — E \ B ~ V ' 
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Fig. 3. — The error in B — V versus the V magni- 
tude for all systems in the survey. 
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Fig. 4. — The error in B — V versus B — V for all 
systems in the survey. 
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Ic)- Using the absorption ratios, 



AW 

A(V) 



given m 

ICardelli. Clavton. fc Mathisl (Il989h and assuming 
a total to selective extinction ratio of Ry =3.1, 
the ratio §^7377 nas a value of 0.647. 

Figure [3] is the color difference diagram for all 
of the systems in the survey. The solid curve 



shows the main se quence from I Walker 
V-I c < 0.95 and lCousinsl (jl978h for V 
0.95. 



(119851) for 
Ic > 

Because Q is reddening-independent, the 



reddening lines in such a diagram are horizon- 
tal and can be used to estimate the reddening 
of a given star. The horizontal line in Figure 
shows the reddening of a main sequence star with 
(V — Ic)o = 0-79 and the vertical ticks indicate 
increasing reddening in steps of 0.1 in E(B — V). 

One concern with this approach to deter- 
mining the reddening is th at many, perhaps all 
(|Pribulla fc Rucinski I I2OO8) , of these objects have 
companions and this can have a large effect on 
the color difference diagram. If, as is usually the 
case, the companions are cool, low-mass objects, 
they will affect the longer wavelength filters more 
than the shorter ones, and make the systems look 
redder than they otherwise would be. However, 
except in most circumstances, we expect the third 
bodies to be of low luminosities, making the effect 
on the system colors fairly small, perhaps a few 
hundredths of a magnitude. Unfortunately, these 
errors are large enough to greatly diminish the use- 
fulness of the color difference diagram for systems 
only slightly affected by interstellar reddening, as 
we expect most of these systems to be due to their 
brightnesses, and thus nearby distances. 

While not as useful as we had hoped for de- 
termining reddening, the color difference diagram 
does distinguish a few systems as being unusual. 
Figure [3] shows a few systems at large horizontal 
distances from the main sequence curve, such as 
VI 044 Her whose E(V — Ic) value from the color 
difference diagram is about 0.9. The question is 
whether that is caused by large interstellar redden- 
ing, or something else, such as a luminous compan- 
ion. Looking at the expected reddening from the 
neutral hydrogen column density provides a clue. 
Predehl fc Schmitt I (jl995l ) used ROSAT obser- 



vations to determine the relationship between X- 
ray absorption and optical extinction. Again as- 
suming Ry =3.1, their results show that E(B — 
V) = 1.8 X 10~ 22 x N HI . To compute the N HI 



value for individual systems, we used the "nh" tool 
in the "heasarc" sub-package of the FTOOLS soft- 
ware available from the NASA Goddard Space- 
flight Center^ Although this approach to estimat- 
ing the reddening doesn't suffer the same prob- 
lems as the color difference method, it does have 
potential problems. The main issue is that the 
Nhi column density is estimated by interpolat- 
ing between the approximately 0.5 degree sam- 
pling points of the Nhi surveys as described in 
the software documentation. Caution must be ex- 
ercised in using this method in areas where there 
are inhomogeneities in the interstellar medium on 
the scale of the interpolation. With that caveat in 
mind, we can compute the reddening and compare 
it to the value from the color difference method. 

In the case of V1044 Her, we find that the value 
of E(B — V) from estimated from the Nhi column 
density is 0.05, while E(B — V) from the color 
difference method is 0.58. In order for the color 
difference value to be correct, the system would 
have to have a (B — V)o value of about 0.34, mak- 
ing it an early F-type star, but its orbital period 
of about 0.24 days precludes a spectral type that 
early. In this case, the reddening estimate based 
on the Nhi column density seems much more rea- 
sonable, but that leaves unanswered the question 
of why the B — V value, implying T e // w 5000if , 
and the V — Ic value, implying T e ff w 4300K, 
are so disparate. The V — Rc and R — Ic values, 
0.69 and 0.60 respectively, are consistent with the 
V — Ic value, implying that the B — V value is too 
blue. The system was measured on three nights 
within a week of one another, and the three B — V 
values were consistent with one another, with val- 
ues of 0.91, 0.92 and 0.92. The images show no 
o bvious problems and t he un filtered light curve 
of lBlattler fc Diethelml (|200ll ) does not appear to 
have large amounts of third light, so the reason for 
the disparity remains unclear. Perhaps the sys- 
tem is more active than usual, leading to elevated 
brightness in the B filter. Being one of only a few 
known systems with an orbital period less than 
0.25 days, further observational study of V1044 
Her is warranted. 

The Nhi method occasionally produces red- 
dening values that are very large compared to the 
color difference method. One example is VI 191 



1 Available at http: / /heasarc .gsfc .nasa. gov/lheasof t/ f tools/heasarc .html 
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Cyg whose E(B — V) from the color difference 
method is 0.05 while the Nhi method yields a 
value of 1.13. This difference is not very surpris- 
ing since Cygnus lies along the galactic plane and 
large inhomogeneities in the interstellar medium 
are to be expected. Because of the weaknesses of 
each method, the safest approach is to trust only 
reddening values where the two methods are in 
good agreement. There are 321 systems that have 
equal values within the error bars, with the error in 
the Njji method assumed to be 0.02 magnitudes. 
However, of those, many have fairly large errors in 
the values from the color difference method, with 
many having errors significantly larger than the 
value itself. To generate a more useful list, we fur- 
ther culled the list by eliminating systems where 
the difference between the two methods was 0.02 
mag or less. This culling leaves 101 systems where 
the reddening values between the two methods are 
consistent. Table 2 shows how the values compare 
for these 101 systems. 

Figure [6] shows the period-color relation for 
the s ystems in Table 2 along with the systems 
from iRucinski I (1997) with periods less than 0.6 
days. The solid line is the short period blue enve- 



lope (SPBE) from lRucinskil (|l998l ) which can be 
likened to the zero-age main sequence for overcon- 
tact binaries in that systems will evolve towards 
redder colors and longer periods. Overcontact sys- 
tems are thus not expected to be bluer than the 
SPBE. The systems in Figure[5]that are bluer than 
the SPBE by more than one standard deviation 
in the B — V value (3091792, 6133549, 3947987, 
2569021, 9649548, 3962973, 8299111, 6225358, 
and V523 Cas) could be systems that are not over- 
contact systems, systems with a blue companion 
or they may be systems where our estimate of 
the reddening is too large. In the case of V523 
Cas (whose overcontact nature appears sol i dly es- 
tablished), S amec. Faulkner, fc Williams I |2004) 
find B — V = 1.03 to 1.09 for the binary and 
(B — V)o = 0.97 ± 0.04 for the primary compo- 
nent, whereas we find B — V = 1.08 ± 0.02 and 
(B - V) = 0.89 ± 0.03 for the binary. Their anal- 
ysis of the O-C digram results in a third body 
mass of about 0.4M Q , which would not explain 
the "blue excess" of the system, and since there 
is good agreement between their reddened B — V 
and ours, the difference may lie in the determina- 
tion of the reddening, or perhaps intrinsic activity 
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Table 1 

The BVR c I c colors of 606 W Ursae Majoris binaries. 



System 


R.A. 


Docl. 


Period 


Num. 


B — V 


V - R c 


Rc-Ic 


V -I c 


v max 






(° J2000) 


(° J2000) 


(days) 


Obs. 














299755 


1.4917 


78.9081 


0.43821290 


3 


0.562 ± 0.016 


0.332 ± 0.003 


0.348 ± 0.004 


0.680 ± 0.005 


11.83 


12.50 


1557554 


1.6542 


55.4561 


0.27249130 


3 


0.874 ± 0.007 


0.510 ± 0.007 


0.432 ± 0.016 


0.942 ± 0.017 


9.33 


9.64 


V1007 Cas 


2.0125 


51.1342 


0.33201060 


3 


0.724 ± 0.027 


0.459 ± 0.004 


0.437 ± 0.010 


0.896 ± 0.011 


12.03 


12.43 


11936075 


3.9833 


6.7457 


0.40117920 


3 


0.519 ± 0.012 


0.308 ± 0.008 


0.332 ± 0.009 


0.640 ± 0.012 


11.30 


11.71 


9107370 


4.0583 


22.1661 


0.42090500 


3 


0.512 ± 0.009 


0.324 ± 0.015 


0.322 ± 0.014 


0.646 ± 0.021 


11.49 


11.82 


CN And 


5.1260 


40.2261 


0.46279111 


4 


0.516 ± 0.012 


0.277 ± 0.011 


0.278 ± 0.014 


0.555 ± 0.018 


9.70 


10.25 


6335770 


6.3292 


32.4657 


0.54349740 


3 


0.377 ± 0.021 


0.231 ± 0.006 


0.236 ± 0.004 


0.467 ± 0.007 


10.25 


10.56 


9093368 


6.7958 


10.0151 


0.52122360 


4 


0.534 ± 0.013 


0.313 ± 0.011 


0.299 ± 0.014 


0.612 ± 0.018 


10.19 


10.61 


14680632 


7.0875 


-14.8881 


0.40264980 


3 


0.577 ± 0.017 


0.351 ± 0.007 


0.351 ± 0.017 


0.702 ± 0.018 


11.56 


11.91 


254036 


7.1167 


78.9619 


0.31741200 


3 


0.586 ± 0.008 


0.356 ± 0.015 


0.350 ± 0.005 


0.706 ± 0.016 


10.56 


10.93 


CL Cet 


7.2667 


-17.2169 


0.62161800 


3 


0.389 ± 0.007 


0.221 ± 0.013 


0.233 ± 0.007 


0.454 ± 0.015 


9.88 


10.00 


3667194 


7.6417 


39.4564 


0.41711880 


3 


0.571 ± 0.018 


0.337 ±0.003 


0.354 ± 0.002 


0.691 ± 0.004 


12.08 


12.47 


9141373 


7.6917 


14.5461 


0.34255000 


4 


0.589 ± 0.011 


0.384 ± 0.017 


0.345 ± 0.017 


0.729 ± 0.024 


11.22 


11.55 


9124077 


8.5542 


20.8739 


0.34014730 


3 


0.596 ± 0.008 


0.392 ± 0.017 


0.368 ± 0.001 


0.760 ± 0.017 


11.91 


12.21 



Note. — Table 1 is published in its entirety in the electronic edition of the Astronomical Journal. A portion is shown here for guidance regarding 
its form and content. 



in the system that enhanced the brightness in the 
B filter during our observations. We note that 
both the color difference and the Njji methods 
show good agreement for the reddening in V523 
Cas (see Table 2). Further observational studies 
of these systems will be needed to determine their 
nature. 

6. Conclusions 

Standardized photometry in Johnson- Cousins 
BVRcIc filters was obtained for 606 known or 
suspec ted W UMa eclipsing binaries. 
( 20061 ). from which a large fraction of the targets 
in this survey were taken, estimate that perhaps 
5% of the systems in their catalog may be ob- 
jects other than W UMa systems, so we have a 
reasonably pure sample of such binaries. Com- 
parison of our results with previously published 
photometry on a few dozen systems shows good 
agreement. Estimates of the interstellar redden- 
ing were made with the color difference method 
and by using the neutral hydrogen column density 
along the line of sight to the stars. While neither 
approach is ideal, they each have strengths and 
are based on independent observational data. A 
subset of 101 systems in the survey shows good 
agreement between the two methods. It is hoped 
that these data will be useful in further research 
to understand the nature of these very interesting 
binaries. 



We thank the referee, Slavek Rucinski, for a 
helpful review that greatly improved the quality 
of the paper. This research has made use of the 
SIMBAD database, operated at CDS, Strasbourg, 
France. 
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Table 2 

Reddening estimates for the 101 systems showing a maximum difference of 0.02 mag in 
e(b - v) between the color difference method and the n hi method. 



System 


E(B - V)cd 


E(B - V)n hi 


System 


E(B - V)cd 


E(B - V)n hi 


System 


E(B - V)cd 


E(B - V) Nhi 


V523 Cas 





.197±0.177 


0.191 


3091792 


0.193±0. 


042 


0.188 


GSC 0067-00348 


0. 


110±0.134 


0.127 


6225358 





.098±0.022 


0.087 


7209961 


0.124±0. 


.042 


0.144 


14541307 


0. 


085±0.031 


0.071 


CW Cmi 





,076±0.042 


0.077 


SX Crv 


0.037±0. 


055 


0.054 


HO Psc 


0. 


054±0.075 


0.054 


3962973 





,118±0.064 


0.106 


5447547 


0.090±0. 


.061 


0.078 


11670840 


0. 


098±0.033 


0.087 


DF Hya 





.030±0.055 


0.047 


V1007 Cas 


0.184±0. 


045 


0.171 


3947987 


0. 


207±0.135 


0.199 


76033 





.079±0.046 


0.081 


13007821 


0.046±0. 


.042 


0.039 


4666562 


0. 


097±0.033 


0.098 


AZ Vir 





,033±0.078 


0.025 


LZ Dra 


0.108±0. 


.082 


0.099 


15697139 


0. 


041±0.053 


0.057 


AC Boo 





.039±0.058 


0.021 


10724466 


0.056±0. 


032 


0.065 


13120541 


0. 


052±0.046 


0.037 


14819714 





064±0.030 


0.080 


14540910 


0.061±0. 


033 


0.071 


CK Boo 


0. 


011±0.033 


0.025 


4701979 





,109±0.051 


0.106 


QW Gem 


0.161±0. 


.064 


0.145 


V829 Her 


0. 


042±0.025 


0.022 


VZ Lib 





,090±0.049 


0.107 


16148635 


0.074±0 


017 


0.085 


5269546 


0. 


008±0.042 


0.017 


V1097 Her 





.052±0.023 


0.059 


AW CrB 


0.019±0 


.080 


0.018 


12209034 


0. 


131±0.056 


0.144 


4750697 





.095±0.049 


0.079 


7329879 


0.052±0 


.040 


0.060 


11899139 


0. 


065±0.024 


0.058 


DZ Psc 





054±0.130 


0.041 


GSC 2545-00970 


0.000±0 


.033 


0.016 


12049671 


0. 


048±0.046 


0.044 


5261236 





.024±0.057 


0.025 


XY Boo 


O.OlliO 


.043 


0.026 


14378336 


0. 


070±0.079 


0.072 


17385764 





.038±0.066 


0.028 


HV Aqr 


0.059±0 


059 


0.062 


7493998 


0. 


037±0.060 


0.025 


1033806 





,055±0.059 


0.041 


LO And 


0.140±0. 


031 


0.155 


11370982 


0. 


110±0.079 


0.130 


7269084 





.086±0.056 


0.093 


BI CVn 


0.017±0 


.021 


0.018 


10142767 


0. 


041±0.030 


0.054 


5214166 





.041±0.036 


0.026 


V1036 Her 


0.069±0 


.061 


0.059 


7336416 


0. 


068±0.092 


0.065 


17357108 





,048±0.021 


0.047 


609692 


0.131±0 


.095 


0.137 


2643685 


0. 


013±0.063 


0.027 


AR CrB 





,043±0.069 


0.050 


908512 


0.070±0 


.046 


0.084 


9221283 


0. 


074±0.033 


0.067 


BL Dra 





.083±0.079 


0.073 


14591442 


0.051±0 


.032 


0.048 


6099330 


0. 


178±0.089 


0.172 


10476945 





.008±0.022 


0.026 


2888729 


0.041±0 


.046 


0.035 


4805949 


0. 


048±0.060 


0.056 


2657877 





024±0.035 


0.021 


3667194 


0.074±0 


.045 


0.075 


15779695 


0. 


079±0.081 


0.091 


10608583 





.053±0.038 


0.055 


4873888 


0.038±0 


.059 


0.020 


4760150 


0. 


072±0.059 


0.059 


2795010 





.023±0.048 


0.023 


13259318 


0.042±0 


.059 


0.028 


13363862 


0. 


032±0.033 


0.050 


XY LMi 





.021±0.060 


0.021 


ET Psc 


0.057±0. 


.040 


0.061 


6375021 


0. 


078±0.046 


0.079 


16134469 





,074±0.094 


0.093 


11905222 


0.041±0 


.023 


0.046 


9311873 


0. 


241±0.041 


0.228 


15595338 





.057±0.142 


0.057 


V728 Her 


0.006±0 


.046 


0.024 


14560232 


0. 


059±0.069 


0.071 


V2759 Ori 





,254±0.047 


0.261 


10141474 


0.054±0 


.042 


0.056 


14304024 


0. 


084±0.104 


0.097 


1206915 





102±0.076 


0.101 


1033721 


0.057±0 


.093 


0.039 


3791537 


0. 


101±0.052 


0.116 


RU UMi 





,037±0.062 


0.019 


V843 Mon 


0.241±0 


.083 


0.224 


BC Eri 


0. 


053±0.130 


0.072 


TZ Lyr 





.043±0.066 


0.052 


ES UMa 


0.056±0 


.079 


0.073 


14678045 


0. 


047±0.052 


0.044 


6335770 





073±0.093 


0.070 


2962451 


0.040±0 


127 


0.043 
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